Quantum Griffiths singularities in ferromagnetic metals 
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We present a theory of the quantum Griffiths phases associated with the ferromagnetic quantum 
phase transition in disordered metals. For Ising spin symmetry, we study the dynamics of a single 
rare region within the variational instanton approach. For Heisenberg symmetry, the dynamics of 
the rare region is studied using a renormalization group approach. In both cases, the rare region 
dynamics is even slower than in the usual quantum Griffiths case because the order parameter 
conservation of an itinerant ferromagnet hampers the relaxation of large magnetic clusters. The 
resulting quantum Griffiths singularities in ferromagnetic metals are stronger than_ power laws. For 
example, the low-energy density of states p(e) takes the asymptotic form exp[{— Alog(eoA)} ]/e 
with A being non-universal. We contrast these results with the antiferromagnetic case in which the 
systems show power-law quantum Griffiths singularities in the vicinity of the quantum critical point. 
We also compare our result with existing experimental data of ferromagnetic alloy Ni^Vi-a,. 

PACS numbers: 71.27.+a. 75.10.Nr, 75.40.-s, 75.50.Cc 
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I. INTRODUCTION 

The low-temperature behavior of quantum many- 
particle systems can be sensitive to impurities, defects or 
other kinds of quenched disorder. This effect is especially 
important near quantum phase transitions, where fluctu- 
ations in time and space become connected. The inter- 
play between static disorder fluctuations and large-scale 
quantum fluctuations leads to much more dramatic ef- 
fects at quantum phase transitions than at classical phase 
transitions, including quantum Griffiths singularities j» — 
infinite-randomness critical points featuring exponential 
rather than power-law scaling^ and the smearing of the 
transition^ 

The Griffiths effects at a magnetic phase transition in 
a disordered system are caused by large spatial regions 
(rare regions) that are devoid of impurities and can show 
local magnetic order even if the bulk system is globally 
in the paramagnetic phase. The order parameter fluctu- 
ations induced by rare regions belong to a class of ex- 
citations known as instantons. Their dynamics is very 
slow because flipping the rare region requires a coherent 
change of the order parameter over a large volume. Grif- 
fiths showed^ that this leads to a singular free energy 
not just at the transition point but in a whole param- 
eter region which is now known as the Griffiths phase. 
In classical systems, the contribution of the rare regions 
to thermodynamic observables is very weak. However, 
due to the perfect disorder correlations in (imaginary) 
time, Griffiths effects at quantum phase transitions are 
enhanced and lead to power-law singularities in thermo- 
dynamic quantities (for reviews see, e.g., Refs. [7| and 
I- 

The systems in which quantum Griffiths behavior was 
originally demonstrated 2 - - — all have undamped dynam- 
ics (a dynamical exponent z = 1 in the clean system). 
However, many systems of experimental importance in- 
volve superconducting 9 or magnetic^ - — degrees of free- 
dom coupled to conduction electrons. This leads to over- 



damped dynamics characterized by a clean dynamical ex- 
ponent z > 1. Studying the effects of the rare regions in 
this case is therefore an important issue. It has been 
shown that metallic Ising antiferromagnets can show 
quantum Griffiths behavior at higher energies, where the 
damping is less important^ In contrast, the quantum 
Griffiths singularities in Heisenberg antiferromagnets are 
caused by the dissipation and occur at lower energies*^ 

In recent years, indications of quantum Griffiths phases 
have been observed in experiments on a number of 
metallic systems suchs as magnetic semiconductors ^r— 
Kondo lattice ferromagnets 19,20 and transition metal 
ferromagnetsi 21 ' 22 All these experimental observations of 
quantum Griffiths phases are in ferromagnets rather than 
in antiferromagnets. However, in contrast to antiferro- 
magnets, a complete theory of quantum Griffiths phases 
in ferromagnetic metals does not yet exist. 

In this paper, we therefore develop the theory of 
quantum Griffiths effects in ferromagnetic metals with 
both Ising and Heisenberg symmetries. We show that 
the quantum Griffiths singularities do not take power- 
law form, in contrast to those in antiferromagnets.^ 
The rare-region density of states behaves as p(e) ~ 
exp[{— Alog(eo/e)} 3 / 5 ]/e in the low-energy limit, where 
A plays a role analogous to the non-universal Griffiths 
exponent. This means, that the Griffiths singularity is 
stronger than a pure power law. This kind of density 
of states leads to non-power-law dependencies on the 
temperature T of various observables, including the spe- 
cific heat, C ~ exp[{-Alog(T /T)} 3 / 5 ], and the mag- 
netic susceptibility, X ~ exp[{-A log(T /T)} 3 / 5 ]/T. The 
zero-temperature magnetization-field curve behaves as 

M ~exp[{-Alog(iViO} 3/5 ]. 

The paper is organized as follows. In Sec. II, we intro- 
duce the model: Landau-Ginzburg- Wilson order param- 
eter field theories for ferromagnetic Ising and Heisenberg 
metals. In Sec. Ill, we study the dynamics of a sin- 
gle rare region. For the Ising case, we use a variational 
instanton calculation, and for Heisenberg symmetry, we 
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use a renormalization group theory of the quantum non- 
linear sigma model with a damping term. In Sec. IV, we 
average over all rare regions and calculate observables 
in the ferromagnetic quantum Griffiths phase. In Sec. 
V, we compare our predictions with existing experimen- 
tal data. Finally, we conclude in Sec. VI by discussing 
the difference between ferromagnetic and antiferromag- 
netic quantum Griffiths singularities as well as some open 
questions. 



stems from the coupling to the conduction electrons. 
In Eq.@, r m is a microscopic time, and in Eq.©, 7 
parametrizes the strength of the dissipation. </>(q, u„) is 
the Fourier transform of the order parameter <f>(r, r) in 
momentum and Matsubara frequency. The value of a 
depends on the character of the electron motion in the 
system and equals I or 2 for ballistic and diffusive ferro- 
magnets, respectively. 



II. THE MODEL 

Rare region effects in disordered metallic systems are 
realized both in Ising magnets^ 4 - and in Heisenberg 
magnets^ In the following, we consider both cases. Our 
starting point is a quantum Landau-Ginzburg- Wilson ac- 
tion of the itinerant ferromagne t 23 ' 24 ' 25 



S — S'stat + S'diss + Sdyn 

where the static part has the form 



(f) 



--E / dr 
Jo 



<0 2 (r,r) + [V^(r,r)] s 



(2) 



Here, Eq is a characteristic energy (assumed to be of the 
order of the band width in a transition metal compound 
or the order of the Kondo-temperature in an /-electron 
system) . We measure lengths in units of the microscopic 
length scale £o- t > is the bare distance of the bulk 
system from criticality. </>(r, t) is the dimensionless order 
parameter field. It is a scalar for the Ising model, while it 
has three components (4>i, 4>2, <j>^) f° r a Heisenberg mag- 
net. 

We consider disorder coupled to the square of the order 
parameter. The corresponding action has the form 



S'd 



--E / dr 
Jo 



d 3 r V(r)0 2 (r,r), 



(3) 



where ^(r) is the disorder potential. 



The dynamical part of the action (p} is Sdyn = S^J n - 



S dyn> where 



S$ a = E T% l I dr I rf a r[9 r 0(r,r)] 



(4) 



corresponds to the undamped dynamics of the system 
with the clean dynamical exponent z = 1, while 



? (2) 

'dyn 



Eo^ 



. I^(q,^n)| 2 



(5) 



describes overdamped dynamics with conserved order pa- 
rameter (clean dynamical exponent z — 2 + a), which 



III. DYNAMICS OF A SINGLE RARE REGION 

In this section, we study the dynamics of a single 
droplet formed on a rare region of linear size L. This 
means, we consider a single spherical defect of radius L 
at the origin with potential V(r) = —V for r < L, and 
V(r) = otherwise. We are interested in the case V > 0, 
i.e., in defects that favor the ordered phase. 

The effective dimensionality of the model defined by 
Eq.([T]) is d c s = 3 + z. Thus, the clean model ([T]) is 
above its upper critical dimension (d c = 4), implying 
that mean-field theory is valid. The mean-field equation 
for a static order parameter configuration 4>o{r) is 26 



V 2 Mr) + [t + V(r)]Mr) + $(r) =0, 



with solution 



Mr) = 



for 
for 



r < L 
r > L. 



(6) 



(7) 



This implies that the order parameter is approximately 
constant in the region r < L and decays outside of it. 

To study the dynamics of the droplet, we start form 
the variational instanton approach. 27 In the simplest case 
the droplet maintains its shape while collapsing and re- 
forming. In order to estimate the action associated with 
this process we make the ansatz 



(r,r) = <f>' (r)r)(T) . 



(8) 



Here, <f)' (r) must be chosen such that J d 3 rtp(r, r) is time 
independent because of order parameter conservation in 
an itinerant ferromagnet. This can be done by introduc- 
ing 0g(r) = 0o( r )(l — Ar) such that the q = Fourier 
component is cancelled. A is a constant to be determined. 
In the limit of a large rare region, Lt » I, we find 



Mr) -Mr) (l-™) 



(9) 



In the following subsections, using ansatz ©, we sep- 
arately discuss the dynamics of the droplet in itinerant 
Ising and Heisenberg ferromagnets. 



A. Itinerant Ising model 

We now calculate the tunneling rate between the "up" 
and "down" states of a single rare region in an itiner- 
ant Ising ferromagnet by carrying out variational instan- 
ton calculations ! 27 ' 28 To estimate the instanton action, 
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we use the ansatz ((8]) (which provides a variational up- 
per bound for the instanton action) with 77 (r) = ±1 for 
t — > ±00. Inserting this ansatz into the action ([lj and in- 
tegrating over the spatial variables yields, up to constant 
prcf actors, 



and 



s. 



(1) 

dyn 



dr[-277 2 (r) + t? 4 (t) 



L 3 / dr[drV(T) 



(10) 



(11) 



The part of the action corresponding to the over- 
damped dynamics becomes 



Q 

>->, 



(2) 
dyn 



dr dr' 



(12) 



where the dimensionless dissipation strength a ~ jL 3+a . 
In order to estimate the action ([TO]) to (fT2|) , we make the 
variational ansatz 



dn _ 26>(t 2 - At 2 ) 
dr t 



(13) 



Summing all contributions, we obtain the instanton ac- 
tion 



S ~ L7to + L't + 1 L 3+a logfo/r, 



(14) 



Minimizing this action over the instanton duration 
gives ro ~ L~ a j^f. Correspondingly, the action is S ~ 
7L 3+a . Then, the bare tunneling rate or tunnel splitting 
behaves as 



e ~ exp(— S) ~ exp(— const, x 7L 3 



(15) 



Thus, the bare tunneling rate decays exponentially with 
iJi+a m ^- ne itinerant Ising ferromagnet unlike the tunnel- 
ing rate in the itinerant Ising antiferromagnet ) 26 i which 
decays exponentially with L 3 . The extra factor L a can 
be understood as follows. To invert the magnetization 
of a rare region of linear size L, magnetization must be 
transported over a distance of order of L, because the 
order parameter conservation prevents local spin flips. 
The rare region dynamics thus involves modes with wave 
vectors of the order of q ~ 1/L. Since the part of the 
action corresponding to the overdamped dynamics is 
inversely proportional to momentum q a , we obtain an 
extra factor L a in the action (fT2"]l . 

Within renormalization group methods, 29 the 
instanton-instanton interaction renormalizes the zero- 
temperature tunneling rate to 



,1/(1-") 



(16) 



This implies that at zero temperature, the smaller rare 
regions with a < 1 continue to tunnel with a strongly 
reduced rate, while the larger rare regions (a > 1) 
stop to tunnel and behave classically, leading to super- 
paramagnetic behavior. 



B. Itinerant Heisenberg Model 

A particularly interesting case are itinerant Heisenberg 
ferromagnets because quantum Griffiths phases have 
been observed experimentally in these systems^ - — We 
now study the dynamics of a single rare region in an itin- 
erant Heisenberg ferromagnet. We make the ansatz 



(r>r) = #,(r)n(T) , 



(17) 



Here, n(r) is a three-component unit vector. After sub- 
stituting (fT7j) into the action (|T|) and integrating over the 
spatial variables, we obtain 



S~gr 2 m dT[dMr)f 



drdr 



, n(r)-n(r') 

(r - t'Y + t, 



(18) 



where the dimensionless coupling constant g ~ L 3 and 
a ~ jL 3+a as before. Because there is no barrier in a 
system with continuous order parameter symmetry, the 
static part of the action is constant. Therefore, we can- 
not solve the problem within the variational instanton 
approach. Instead, rotational fluctuations must be taken 
into account. 

We calculate the characteristic relaxation time of the 
rare region by a renormalization group analysis of the 
action (j 18[) . As shown in the Appendix, for weak damp- 
ing a < j, there are two different regimes, where the 
behaviors of the relaxation times are different. Partic- 
ularly, for energies to larger than some crossover energy 
ijj c ~ a/g, undamped dynamics is dominant, and the re- 
laxation time of the rare region has the form 



So 



(19) 



which leads to a power-law dependence of the rare-region 
characteristic energy on L, 

e~L-' 3 . (20) 

For energies uj <S uj Ci overdamped dynamics dominates 
the system properties, and the relaxation time of the rare 
region behaves as 



C ~ exp [const, x ^L i+a ] . 
This results in a characteristic energy of 
e ~ cxp[— const, x 7-L 3+a ] . 



(21) 



(22) 



Thus, the behavior of the characteristic energy in the 
itinerant Heisenberg magnet is the analogous to that of 
the tunneling rate in the Ising model discussed above. 

We can now roughly estimate the size L c of the rare 
region corresponding to the crossover of the two regimes. 
By comparing Ea. (TT9l and Eg . (12T1) . we find for small a: 



(23) 



[log(const./ 7 )/7] 1/(3+Q) . 
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For small rare regions, L < L c , the undamped dynamics 
dominates systems properties and the characteristic en- 
ergy is given by Eq . (f20| . while for L > L c , the damping 
term is dominant and the characteristic energy is deter- 
mined by Ed.([22]>, 

For large damping a 3> <?, the overdamped dynamics 
dominates the system properties for all energies us. Cor- 
respondingly, the characteristic energy is given by (|22|) . 



IV. OBSERVABLES 

In the last section, we have seen that metallic Ising fer- 
romagnets display modified Griffiths behavior at higher 
energies (fT5|) . while at asymptomatically low energies, the 
rare regions freeze and lead to a smeared phase transition 
(|16[) . For Heisenberg ferromagnets, we have found con- 
ventional behavior at higher energies (l20l) . and modified 
Griffiths behavior at low energies (|22[) . Correspondingly, 
we expect modified Griffiths singularities in thermody- 
namic quantities at low energies for itinerant Heisenberg 
ferromagnets, while for metallic Ising ferromagnets they 
should occur at higher energies. 

In this section, we use the single-rare-region results of 
Sec. Ill to study the thermodynamics in these ferro- 
magnetic quantum Griffiths phases. To do so, we need 
to estimate the rare-region density of states. By basic 
combinatorics (see, e.g., Refs. andU), the probability 
for finding an impurity-free rare region of volume L 3 is 
V ~ exp(— bL 3 ) with b being a constant that depends on 
the disorder strength. Combining this and Eq.(f22"j) gives 
the density of states (of the Heisenberg system) in the 
low-energy regime as 



p(e) ~ -exp[{-Alog(e /e)} 



3/(a+3)i 



(24) 



Here, eo is a microscopic energy scale, and the non- 
universal exponent A ~ 5( a + 3 )/3^ pl avs a role similar 
to the usual quantum Griffiths exponent. The same den- 
sity of states follows from Eq. (fT"5l) for the higher-energy 
regime of the Ising model. Thus, in ferromagnetic metals, 
the rare-region density of states does not take power-law 
form, in contrast to the one in antiferromagnets. 

We can now find observables using the rare-region den- 
sity of states (|2~41 . The number n of free rare regions at 
temperature T behaves as 

n(T)~ J dep{e)e- e ' T /{l + e" e / T ) 

~exp[{-Alog(T /T)} 3 /( a + 3 )], (25) 

where T is a microscopic temperature scale. 

The uniform static susceptibility can be estimated by 
summing Curie susceptibilities for all free rare regions, 
yielding 

X (T) = n(T)/T^iexp[{-Alog(T /T)} 3 /( a + 3 )] .(26) 



The dependence of the moment p of the rare region on 
its energy leads to a subleading correction only. 

The contribution of the rare regions to the specific heat 
C can be obtained from 



AE 



dep{e)e e~ e/T /(l + e" e/J ) 



Texp[{-Alog(T /T)} 3 /( Q + 3 )] , 



(27) 



which gives AC* ~ exp[{-Alog(T /T)} 3 /( a+3 )]. Knowing 
the specific heat, we can find the rare region contribution 
to the entropy as AS ~ exp[{-Alog(T /T)} 3 /( Q+3 )]. 

To determine the zero-temperature magnetization in 
a small ordering field H , we note that rare regions with 
e < H are (almost) fully polarized while the rare regions 
with e > H have very small magnetization. Thus, 



dep(e) ~ exp[{-Alog(#o/#)} 



3/(o+3)i 



(28) 



where Hq is a microscopic field (again, the moment of the 
rare region leads to a subleading correction). The zero- 
temperature dynamical susceptibility can be obtained by 
summing the susceptibilities of the individual rare regions 
using the density of states ([2"4")l . 



X(w) = / de/?(e)xrr(w;e) 
Jo 



(29) 



where the dynamical susceptibility of a single rare region 
in Heisenberg metals at zero temperature is given by 30 



Xrr(w + iO;e) 



IJLU 



(30) 



where p is the moment of the rare region. Substituting 
Eq.p)]) into (J29J) we find 



x(c + z0)~ (1+ 7f (a;)) ex P [{-AlogN/a;|} 3/( " +3) ] ; 

(31) 

where luq is a microscopic frequency. This result can be 
used to estimate the rare region contribution to the NMR 
spin relaxation time T\. Inserting pip into Moriya's 
formula 31 for the relaxation rate yields 



T 



l/Ti ~ ^exp[{-Alog|cVw|} 3/(a+3) ] 



(32) 



EXPERIMENT 



Recently, indications of a quantum Griffiths phase 
have been observed in the transition metal ferromagnet 
Ni a Vi-3;j 21 i 22 The behavior of the thermodynamics has 
been described well in terms of the power-low quantum 
Griffiths singularities predicted for an itinerant antifcrro- 
magnet (and the transverse-field Ising model). Here, we 
compare our new theory of ferromagnetic quantum Grif- 
fiths phases with the experimental data given in Refs. 
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FIG. 1. (Color online). Temperature dependence of the sus- 
ceptibility of Ni^Vi-i, for different Vanadium concentrations. 
Solid and dotted lines represent fits to the Eg. (|26p in the 
different temperature r ang es 5 to 300K and I to 20K, respec- 
tively, (data from Ref. \2w 

l2ll and I22I. The residual resistivity of Ni^Vi-^ close to 
the quantum phase transition is rather flight Thus, we 
choose a = 2 for a diffusive ferromagnet. Fig. Q] shows 
the behavior of the susceptibility as a function of tem- 
perature. The curves corresponding to the concentra- 
tions x = 13.0% and x = 15.0% (which are far away 
from the critical concentration x c « 11.5%) are described 
better by power laws rather than our modified quantum 
Griffiths behavior (|26|) . at least above T w 10K (the low- 
temperature upturn is likely due to freezing of the rare re- 
gions). For concentrations x = 12.07% and x — 12.25%, 
our theory fits better than power-law Griffiths singular- 
ities and extends the fit range from 30-300K down to 
5-300K. The curves corresponding to the concentrations 
x = 11.4% and x = 11.6% can be fitted by Griffiths 
power-laws only in the temperature range 30 to 300K, 
our new functional form (126|) does not improve the fit of 
these curves. 

We also compared the prediction (|2"51) for a modified 
magnetization-field curve with the data given in Refs. [2l| 
and [U We found that the fits to power-laws and to 
the modified quantum Griffiths behavior |28|) cannot be 
distinguished. 

Let us also point out that the susceptibility data in the 
temperature range below 20K can be fitted reasonable 
well by Eq. (f26]) . see details in Fig. [TJ Further exper- 
iments may be necessary to decide whether our theory 
applies in this region. 

Overall, our theory does not significantly improve the 
description of the data of Refs. Hi] and over the tem- 
perature range where Griffiths behavior is observed. A 



possible reason is that the relevant rare regions are too 
small. At concentrations x — 13.0% and x = 15.0%, they 
have moments of about fi w 5/ub and // « 1/^b, respec- 
tively. Correspondingly, the effect of the order parameter 
transport cannot play any role, whereas our functional 
forms arise for large rare regions where the order param- 
eter transport limits the relaxation of the rare region. 
A possible reason why the curves corresponding to the 
concentrations x = 11.4% and x — 11.6% can not be 
described by our theory at T < 30K might be that the 
curves are actually slightly on the ordered side of the 
quantum phase transition. 



VI. CONCLUSIONS 

In summary, we studied the dynamics of rare regions 
in disordered metals close to the ferromagnetic quantum 
phase transition, considering the cases of both Ising and 
Heisenberg spin symmetries. The overall phenomenology 
is similar to the well-studied antiferromagnetic quantum 
Griffiths behavior . 6 i 14 ' 15 i 27 Namely, for Ising symmetry 
at low temperatures, the overdamping causes sufficiently 
large rare regions to stop tunneling. Instead, they behave 
classically, leading to super-paramagnetic behavior and a 
smeared quantum phase transition. In contrast, at higher 
temperatures but below a microscopic cutoff scale, the 
damping is unimportant and quantum Griffiths singular- 
ities can be observed. In contrast to the Ising case, the 
itinerant Heisenberg ferromagnet displays quantum Grif- 
fiths singularities when damping is sufficiently strong, 
i.e., at low temperatures. Above a crossover tempera- 
ture, conventional behavior is expected. 

Although the phenomenologies of the ferro- and an- 
tiferromganetic cases are similar, the functional forms 
of the quantum Griffiths singularities are different. In 
ferromagnetic quantum Griffiths phases, the tunneling 
rate (or characteristic energy) of a rare region decays as 
exp[— const, x jL a+3 ] with its linear size L, where a is 
equal to 1 or 2 for ballistic and diffusive ferromagnets, 
respectively. This leads to the modified non-power law 
quantum Griffiths singularities in thermodynamic quan- 
tities, discussed in Sec. IV, in contrast to the power-law 
quantum Griffiths singularities in itinerant antiferromag- 
nets. The reason is the following. Because of the order 
parameter conservation in the itinerant quantum ferro- 
magnet, the damping effects are further enhanced as the 
dimensionless dissipation strength a for a rare region of 
linear size L is proportional to L a+3 rather than L 3 . 

In strongly disordered system, where our theory is 
most likely to apply, the motion of the electron is diffu- 
sive. Correspondingly, we expect a = 2. In hypothetical 
systems with rare regions, but ballistic dynamics of the 
electrons, a would take the value 1. 

In our explicit calculations, we have used Hertz's 
forrn^ of the order-parameter field theory of the itin- 
erant ferromagnetic quantum phase transition. However, 
mode-coupling effects in the Fermi liquid lead to an ef- 
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fective long-range spatial interaction between the order 
parameter fluctuations^ - — In the order-parameter field 
theory, this leads to a nonanalytic momentum depen- 
dence of the static action @. The effects of this long- 
range interaction on the existence and energetics of a 
locally ordered rare region were studied in detail in Ref. 
[2H . This work showed that the long-range interactions 
only produce subleading corrections to the droplet free 
energy. Therefore, including these long-range interac- 
tions in the action ([T]) will not change the results of the 
present paper. 

Let us now turn to the limitations of our theory. In 
our calculations, we assumed that the droplet maintains 
its shape while collapsing and reforming. Correspond- 
ingly, our calculation provides a variational upper bound 
for the instanton action. There could be faster relax- 
ation processes, however, it is hard to image the droplet 
dynamics to avoid the restriction coming from the or- 
der parameter conservation. We treated the individ- 
ual, locally ordered rare regions as independent. But, 
in a real metallic magnet, they are weakly coupled by 
a Rudcrman-Kittel-Kasuya-Yosida (RKKY) interaction 
which is not included in the Landau- Ginzburg- Wilson 
action (TIJ. At the lowest temperatures, this RKKY in- 
teractions between the rare regions induce a cluster glass 
phased Finally, our theory does not take the feedback 
of the order parameter fluctuations on the fermions into 
account. It has been found that for some quantum phase 



transitions, the Landau-Ginzburg- Wilson theory breaks 
down sufficiently close to the transition point due to this 
feedback. 37 ' 38 For strongly disordered systems, this ques- 
tion has not been addressed yet, it remains a task for the 
future. 

Turning to experiment, our theory does not sig- 
nificantly improve the description of the data of 
Ni x Vi_a;j 21 i 22 We believe that the main reason is that our 
theory is valid for asymptomatically large rare regions 
where the order parameter transport plays an important 
role, whereas the experimental accessible rare regions in 
Ni^Vi-a; are not large enough for the order parameter 
conservation to dominate their dynamics. We expect our 
theory can be applied in systems where one can observe 
Griffiths singularities at lower temperatures leading to 
larger rare regions. 

Appendix: Renormalization group theory 

In this Appendix, we show the derivation of 
Eqs. (|I9|21|) by renormalization group (RG) analysis. At 
low temperatures, the action (|18j) is formally equiva- 
lent to a quantum non-linear sigma model 39 in imagi- 
nary time r. We can set n(r) = (it(t),<t(t)), where 
tt(t) = (tt\(t), 7T2(t)) represents transverse fluctuations. 
After expanding in ir and keeping terms up to 0(g~ 2 ), 
0(a- 2 ), we find 3 ^ 



S = 



j ^7T i 9U}2 + I !^ W )' 2 + / ~~72^~^ (f ^ ~ gUJlUJ3 ) np( u l)npMni3'(U3)ni3>(-Ul -U2- ^3) ■ 

(A.l) 



We now consider the case of the small damping aCg. 
Two different energy regimes can be distinguished: (i) 
u> larger than some crossover energy ui c ~ a/g, imply- 
ing that the undamped dynamics dominates the systems 
properties, and (ii) u> <C uj Ci when the damping term is 
dominant. 

(i) Because the contribution of the undamped dynam- 
ics is dominant in this regime, we neglect the damp- 
ing term and renormalize g. To construct a pertuba- 
tive renormalizaition group, consider a frequency region 
[—A, A] (A is a high energy cut off), and divide the modes 
into slow and fast ones, tt(uj) = tt < (uj) + 7r > (w). The 
modes tt < (uj) involve frequency —A/b < uj < A/6, and 
are kept. We integrate out the short- wavelength fluctu- 
ations 7r > (o;) (with frequencies in the region —A < uj < 
—A/b and A/b < uj < A) in perturbation theory using 
the propagator (tt^ (uj)tt^, (uj 1 )) = n5pp'6(co + u')/(guj 2 ). 

After applying standard techniques, we find that this 
coarse graining changes the coupling constant g to g co = 
g + I g (b), where I g (b) = (27rA) _1 (6 - 1). After rescaling 
t' = r/b and renormalizing 7t'(t') = 7r < (T)/C g , we obtain 



the renormalized coupling constant in the form 

a' = b-\ 2 g g co ■ (A.2) 

To find the rescaling factor £ s , we average n over the 
short wavelength modes n > and obtain 

<n)> =<(7T 1 < + 7T> V /l-(^ < +^ > ) 2 )) > 

=(1 - <(tt>) 2 )>/2 + 0(g- 2 ))(n<, v^W) . 

(A.3) 

Thus, we identify 

C a = 1 - <(7r>) 2 >>/2 + 0(g- 2 ) = 1 - ^ + 0(g- 2 ) . 

(A.4) 

Correspondingly, the renormalized coupling constant 
given in Eq. (|A.2[) becomes 

g' = b-\g-I g (b)). (A.5) 
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Setting b — 1 + SI, and integrating Eq. (|A.5[) gives the 
recursion relation g(l) — g(0)e~ l . To find the relaxation 
time, we run the RG to g(l) — 1 and use £ r ~ e l . This 
gives 

C g ~ L 3 . (A.6) 

(ii) In the same way, for low energies w <C uj c , we ne- 
glect the term corresponding to the undamped dynamics 
and renormalize the a coefficient. We find that a is not 
modified by the perturbation, i.e. a co — a and the field 
rescaling factor £ a is given by 

Ca = l_Z^ +0(a -2 )) (A . 7) 



where I a (b) — 27r~ 1 log(&). Then, we find the recursion 
relation a(l) — a(0) —Air^ 1 !. This leads to the relaxation 
time 

a ~ exp[const. x 7 L 3+Q ] . (A.8) 
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